Multiband tight-binding model of local magnetism in Gai^Mn^As 
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We present the spin and orbitally resolved local density of states (LDOS) for a single Mn impurity 
and for two nearby Mn impurities in GaAs. The GaAs host is described by a sp 3 tight-binding 
Hamiltonian, and the Mn impurity is described by a local p-d hybridization and on-site potential. 
Local spin-polarized resonances within the valence bands significantly enhance the LDOS near the 
band edge. For two nearby parallel Mn moments the acceptor states hybridize and split in energy. 
Thus scanning tunneling spectroscopy can directly measure the Mn-Mn interaction as a function of 
distance. 

PACS numbers: 75.50.Pp, 75.30.Hx, 75.30.Gw, 71.55.-i 



The successful growth of ferromagnetic diluted mag- 
netic semiconductors (DMSs) based on III-V com- 
pounds 0,00,0,0] has attracted much attention due to 
the nearby metal-insulator transition and also for their 
potential application to nano-scale nonvolatile storage 
devices and to quantum computation 0, 0- I n order 
to turn a nonmagnetic semiconductor into a magnet, a 
sizable amount (~ 1%) of magnetic dopants are intro- 
duced. Experimental studies 

B M III 113 have 
shown strong correlations between the ferromagnetism 
and the hole carriers that are also contributed by the 
Mn doping in the III-V DMSs. The emerging picture of 
this ferromagnetism is that the Mn magnetic moments 
are localized, and the ferromagnetic coupling is medi- 
ated through the delocalized hole carriers jy, 0. A 
deeper understanding of this ferromagnetic coupling has 
been hampered by insufficient knowledge of the valence 
band structure. The assumption that the holes reside 
in the unperturbed valence bands of the host semicon- 
ductors 0, has recently been questioned by both 



theoretical and experimental studies [lfj, [13, lla, UM ■ Ef- 
fective mass theories provide good spectral resolution at 
the band edge, but cannot describe distortions on dis- 
tances of one or two lattice constants for they are only 
accurate for momenta close to the T point. Density func- 
tional calculations [ItI Hif , however, describe local prop- 
erties well, but state of the art supercell calculations have 
not had sufficient spectral resolution to resolve the shal- 
low bound states in the gap and the sharp resonances in 
the valence bands. 

In this Letter we study the local density of states 
(LDOS) of very dilute concentrations of Mn in GaAs as a 
prototype III-V DMS, and show that the valence bands 
are strongly altered by the Mn dopants, and in return 
influence the interaction between the Mn magnetic mo- 
ments. In order to obtain both sufficient spectral resolu- 
tion and proper dispersion relations throughout the full 
Brillouin zone, a multiband tight-binding approach incor- 
porating spin-orbit interaction is employed. Our results 
show that the hybridization between the Mn 3<i-orbitals 
and the GaAs valence bands leads to spin-polarized reso- 
nances within the valence bands and to delocalized ferro- 



magnetic interaction. The LDOS near the valence band 
edge is significantly enhanced by these resonances. Each 
Mn dopant can enhance the LDOS by as much as a factor 
of 2 up to the second- nearest neighbors (See Fig.QJ, cor- 
responding to a ~ 7% increase of the average LDOS for 
1% Mn concentration. The strong LDOS enhancement is 
consistent with recent angle-resolved photoemission mea- 
surements , and qualitatively explains the increases of 




121 as well 
13. In 



the absorption coefficients in the intraband 
as the interband absorption spectroscopies 
particular, the experimental result that the band-edge 
absorption coefficient increases in finite magnetic fields 
more for one polarization of light is consistent with our 
result that the resonances are both orbitally- and spin- 
polarized. 

Whereas the spin-polarized band-edge resonances will 
alter the Mn-Mn interaction, the acceptor states provide 
a direct means of measuring that interaction. The accep- 
tor level of Mn splits when two Mn magnetic moments 
are in a parallel configuration [lij . The range of interac- 
tion between the Mn dopants can be probed by the size 
of the splitting as a function of the separation. We find 
that the splittings are as large as tens of meV even when 
two Mn dopants are separated by a few lattice constants. 
This magnetic interaction is anisotropic with respect to 
the axis connecting the two Mn dopants due to spin-orbit 
interaction, in qualitative agreement with the continuum 
model [io]. Both the local enhancement of the valence 
band edge and the splitting of the acceptor level could be 
probed by scanning tunneling spectroscopy (STS) |2l|. 

As pointed out by Vogl and Baranowski |22T| . the 
fivefold degenerate Mn 3cZ-orbitals split in a cubic lat- 
tice to an E'-symmetric (d x 2_ a 2-like) doublet, which is 
only weakly coupled to the tetrahedral host, and a T2- 
symmetric (d xy -lihe) triplet, which can effectively couple 
to the neighboring dangling sp 3 -hybrids. The antibond- 
ing states of the T2-symmetric states and the sp 3 -hybrids 
form the acceptor states within the gap. These anti- 
bonding states are delocalized in space since they overlap 
strongly with the host valence bands. In our calculations 
the hybridization of the majority Mn d-orbitals is treated 
as an effective extended spin-dependent potential U\ act- 
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ing on the four nearest-neighbor sites of Mn, because only 
s and p orbitals are explicitly used in our tight-binding 
Hamiltonian. The operator form of our Mn potential is 
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where Uq is the on-site orbital energy difference, c] s (R) 
(q s (R)) is the creation (annihilation) operator of a spin- 
s electron in the I orbital at site R. The Mn dopant 
is located at Ro, and the four nearest-neighbor sites are 
labeled by R1-R4. We neglect the difference of the po- 
tential matrix elements between the s and p orbitals, be- 
cause the s-orbitals make a negligible contribution to the 
LDOS near the valence band edge. The quantization axis 
for the spin is aligned with the Mn core spin. The mi- 
nority Mn d-orbitals are assumed to be much higher in 
energy, and the hybridization energy of them with the 
spin-down states is neglected. We find that the effective 
potential U± is of critical importance in inducing the ac- 
ceptor level at the experimentally observed energy. The 
significance of the hybridization to the acceptor level en- 
ergy has also been noted recently by Dietl et at. . As 
for the delocalized nature of the acceptor state, our re- 
sults show that only about 10% of the spectral weight of 
the acceptor state is concentrated at the Mn site, 20% 
is distributed over the four nearest-neighbor sites, and 
the remaining 70% is extended to farther sites. Further- 
more, the large spin-orbit interaction of GaAs splits the 
gap states into three different energy levels, and intro- 
duces anisotropy to the spatial structure. 

We use the Koster-Slater technique [24| to calculate the 
Green's function whose imaginary part gives the LDOS. 
This method has been proven to give the correct chemical 
trend of impurity levels in semiconductors |2*3 |. and we 
have previously applied it to successfully predict STM 
spectra near impurities in superconductors |2fij . Start- 
ing with the tight-binding Hamiltonian Hq (k) of homo- 
geneous GaAs, one first calculates the retarded Green's 
function, C?o(k, uS) = [u> — Ho(k) + iS}^ 1 . To obtain 
a good description of the valence band structure, we 
use the sp 3 model including spin-orbit interaction |27j |. 
Then, by Fourier transforming Go(k, oj), we construct 
the homogeneous Green's function in coordinate space, 
Go(Ri, Rj, uS), where Rj and Rj label the zincblende lat- 
tice sites. This step consumes the majority of computa- 
tion time. It takes about one day for one link (Rj — Rj) 
with a spectral range of 2 eV on a personal computer. 
The number of links used in this study is about 500. A 
constant linewidth of S = 10 meV produces a good spec- 
tral resolution within a reasonable computation time. 
The final Green's function is obtained by solving the 
Dyson's equation, 

g{lo) = GoH + [i-GoHy] _1 GoH , (2) 
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FIG. 1: LDOS spectra for the charged (ionized) Mn dopant in 
Gai-zMn^As (corresponding to bulk n-doped) at (a) the Mn 
site, at (b) one of the nearest-neighbor As sites, and at (c and 
d) one of the second-nearest-neighbor Ga sites. The shaded 
balls in each panel show the corresponding atomic sites, and 
the black ball shows the Mn site. The Mn core spin is aligned 
with the z-axis of the lattice, and the z-axis is aligned ver- 
tically. In this particular case, the spectra at four nearest- 
neighbor sites are identical, and have two different functional 
forms at the second-nearest-neighbor sites. The total LDOS 
is shown by the solid lines. The dashed lines show the projec- 
tions onto the two spin polarizations relative to the Mn core 
spin. The dotted lines show the LDOS of the bulk GaAs. 
The zero of energy is the valence maximum of bulk GaAs. 
The valence band LDOS of bulk GaAs is more concentrated 
at the As sites than at the Ga sites, because the LDOS is 
the density of states modulated by the probability density of 
electrons. For the neutral Mn dopant (corresponding to bulk 
p-doped Gai-^Mn^As) only the highest energy peak of the 
three localized peaks apparent above in each panel would be 
visible. 



where G(u>) is the full matrix representation using all 
atomic orbitals at all lattice sites. The LDOS at each 
site Rj is given by 



A(Ri,uj) = 1m tr a Gr(Rj,Rj,w) 

7T 



(3) 



where the tr a is taken with respect to the orbitals of the 
a atom, depending on which type of atom is actually 
located at the site Rj. 

To determine the values of Uo and U\, we first consider 
the possibility that the potential is nonzero only at the 
Mn site (U\ = 0), and the size of the on-site potential is 
assumed to be the same for all orbitals. When the Mn 
atom replaces a Ga atom, the strength of this on-site po- 
tential is estimated to be Uo ~ 1 eV based on the energy 
difference of the ionization energies of Ga (4s 2 4p) and of 
Mn (3c£ 5 4s4p), which are about 6 and 5 eV respectively. 
Not only is this on-site potential too weak to bind any ac- 
ceptor level in the gap, it is not possible (even with an un- 
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FIG. 2: (color) The spatial structure of (a) the Mn acceptor 
level at 113 meV and of (b) the valence band LDOS at —500 
meV. The LDOS peaks are centered at the Mn atom, and are 
extended mainly on the (001) plane. In (b), far away from 
the Mn site, the sites with higher LDOS are occupied by the 
As atom, and with lower LDOS are occupied by the Ga atom. 
We have assumed that the squared modulus of the Wannier 
functions to be a Gaussian with a width of half the distance 
between neighboring Ga and As atoms. 
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FIG. 3: The orbitally resolved LDOS for (a) p x + ip y , (b) p z 
and (c) p x — ip y orbitals. The Mn core spin is aligned with 
the z-axis. Each panel contains four LDOS curves. From 
top to bottom, the curves show the LDOS at (i) the Mn, (ii) 
the nearest-neighbor (As) sites, and (iii and iv) the second- 
nearest-neighbor (Ga) sites. Cases (iii) and (iv) correspond 
to the configurations (c) and (d) in Fig. Q respectively. The 
LDOS spectra are shifted by multiples of one unit for easier 
visualization. 



realistically large on-site potential, Uq ^> 10 eV) to bind a 
hole more than 60 meV above the valence band edge. The 
Mn dopant introduces an acceptor level at 113 meV above 
the valence band edge, therefore, the effective potential 
must be extended at least to the four nearest-neighbor As 
sites. If we fix Uq = 1 eV, and tune U\ to give the cor- 
rect acceptor level energy, the required nearest-neighbor 
potential is about 3.59 eV. This nearest-neighbor poten- 
tial is too large to be accounted for by the tail of the 
screened Coulomb potential. The screened Coulomb po- 
tential at the nearest-neighbor sites of a singly charged 
Mn center is only about 0.5 eV, or even smaller if the 
system is doped with carriers, as in our case. Therefore, 
the large effective potential at the nearest-neighbor sites 
must result from the hybridization of the Mn 3<i-orbitals 
and the host sp 3 -hybrids. The exchange interaction ef- 
fectively pushes the valence-band states with the same 
spin polarization as the Mn core spin into the gap. The 
states with the opposite spin polarization essentially re- 
main unperturbed. 

Figurenshows the LDOS spectra at the Mn site and at 
neighboring sites. In this particular case we have chosen 
the Mn core spin to align with one of the (100) crystal 
axes (which is the bulk easy axis p^L As a result, the 
spectra at the four nearest-neighbor As sites are identi- 
cal due to the residual symmetry operations of Ci and 
S4 about the crystal axes. The twelve second-nearest- 
neighbor Ga sites are divided into two classes with four 
and eight sites each. The acceptor states are almost fully 



spin-polarized (parallel to the Mn core spin) and split 
into three energy levels due to the spin-orbit interaction. 
For bulk n-doped Gai_ x Mn x As the Fermi level lies above 
the upper acceptor state, so the impurity is ionized and 
all three levels would be visible in a tunneling experi- 
ment. If the Fermi level lies below the upper acceptor 
state, such as in bulk p-doped Gai-^Mn^As, then due to 
the electron-electron interaction only one localized state 
will be visible above the Fermi level and none below. 
Due to band-bending at the surface it is often possible 
to see both the ionized and neutral dopants in the same 
sample (as in Ref. Elj) . As the Mn core spin rotates, 
the orbital character of the acceptor states changes, but 
the three energy levels remain the same within the sp 3 
model. There is some experimental evidence for the ex- 
istence of the two additional acceptor energy levels |2^. 
The significant enhancement of the spin-up LDOS contin- 
ues into the valence band to 1.5 eV below the band edge. 
The spin-polarized resonances emerge at an energy close 
to the split-off band top (<~ 350 meV below the band 
edge). On the other hand, the spin-down band is weakly 
perturbed and non-resonant. Figure H ( c ) an d (d) also 
show that the large enhancement of the valence band 
LDOS extends to the second-nearest neighbors. From 
the spectral weight of the acceptor state, ~ 5% at each 
As site, we estimate the p-d exchange interaction (3Nq to 
be \\J X x 20% « -0.3 eV, within a factor of 3 of that 
obtained from transport measurements (see Chapter 1 in 
Ref. |6|). The discrepancy may come from the enhanced 
LDOS near the Mn spin, which is not included in the 
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FIG. 4: Splitting of the acceptor levels. The left panel shows 
the spectra at the Mn site with two Mn dopants separated 
by (a) (0,0,3) and (b) (-3/2,3/2,0). The dashed line shows 
the spectrum when there is only one Mn dopant. In (b) the 
two split peaks of the 113 meV level are indicated by the two 
heavy-line segments. The right panel shows the amount of 
the splitting of the 113 meV level. One of the Mn atoms is 
located at the (0, 0, 0) site (bottom-left corner), and the other 
is at one of the circles. The number in the circle shows the 
splitting energy in meV. Both Mn core spins are aligned with 
the 2-axis. The black dots show the As sites. 



orbital character more extended in the plane perpendic- 
ular to the Mn core spin, the splitting of the acceptor 
level is the largest when the Mn core spins are oriented 
perpendicular to the axis that joins them. 

In summary we have presented calculations of the 
LDOS near Mn in GaAs. The p-d hybridization in- 
duces spin-polarized resonances that enhance the LDOS 
near the valence band edge. The spin-orbit interaction 
splits the acceptor level, and introduces anisotropy to the 
Mn-Mn interaction in addition to the crystal anisotropy. 
The range of the Mn-Mn interaction is shown to extend 
through several lattice constants. For quantitative com- 
parisons with STS measurements, significant surface ef- 
fects are expected for a Mn atom at the top layer because 
the Mn atom has one less As neighbor, and the hybridiza- 
tion changes. However, it has been reported that do pan t 
atoms as deep as in the fifth layer are seen by STS [2l|. 
and the spectra for these cases should closely resemble 
the bulk. 

This work was supported under the ARO MURI 
DAAD19-01-1-0541. 



transport modeling. 

The spatial structure of the LDOS is further illustrated 
in Fig. [21 For the Mn potential used here, the upper ac- 
ceptor state is peaked at the Mn site. The ratio between 
the peak heights at the Mn and at the nearest-neighbor 
As sites is related to the strength of Uo and U\. Gen- 
erally, increasing the potential strength at one site re- 
duces the spectral weight at the same site. The degree 
of acceptor-state localization apparent in Figgis deter- 
mined largely by the effective band masses of the GaAs 
host, which are dependent on the tight-binding parame- 
ters of Ref. |U If the spin is aligned with the z axis, the 
uppermost acceptor state is spatially extended in the x-y 
plane, because the orbital character, shown in Fig. |3 is 
Px + iPy 

We also calculate the LDOS for two nearby Mn dopants 
by inverting the Dyson equation for a two-impurity po- 
tential. When the two Mn core spins are parallel, the ac- 
ceptor states centered at the two Mn sites interfere, and 
further split into "bonding" and "antibonding" states. 
When the two Mn core spins are antiparallel, one of the 
U\ terms acts on spin-down As-orbitals, and the poten- 
tial for the spin-up (down) states at one Mn site is the 
same as for the spin-down (up) states at the other Mn 
site. Therefore, the acceptor states remain in three de- 
generate levels. The splitting of the acceptor level for 
two parallel Mn core spins can be used as a probe for 
the range of interaction between Mn dopants. Figure 0] 
shows the long-ranged nature of the interaction. It also 
shows the anisotropy of the interaction resulting from the 
spin-orbit interaction. Because the acceptor state has an 
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